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A new generation of upcoming space-based experiments will soon start to probe the spectrum of 
cosmic ray antiparticles with an unprecedented accuracy and, in particular, will open up a window 
to energies much higher than those accessible so far. It is thus timely to carefully investigate 
the expected antiparticle fluxes at high energies. Here, we perform such an analysis for the case of 
antiprotons. We consider both standard sources as the collision of other cosmic rays with interstellar 
matter, as well as exotic contributions from dark matter annihilations in the galactic halo. Up to 
energies well above 100 GeV, we find that the background flux in antiprotons is almost uniquely 
determined by the existing low-energy data on various cosmic ray species; for even higher energies, 
however, the uncertainties in the parameters of the underlying propagation model eventually become 
significant. We also show that if the dark matter is composed of particles with masses at the TeV 
scale, which is naturally expected in extra-dimensional models as well as in certain parameter regions 
of supersymmetric models, the annihilation flux can become comparable to - or even dominate - 
the antiproton background at the high energies considered here. 

PACS numbers: 98.70.Sa, 96.50.sb, 95.35. +d 



I. INTRODUCTION 

Since the pioneering works of Victor Hess almost 100 
years ago we know that the earth is constantly exposed 
to a bombardment of cosmic rays with presumably galac- 
tic origin. Extending up to extremely high energies, they 
mainly consist of protons (90%) and helium (9%), but 
also contain sizeable contributions of other nuclei and 
antiparticles. Most of these particles have probably been 
accelerated in the blast waves of supernova remnants, 
while others - like antiprotons - are rather produced by 
the spallation of the interstellar medium by primary cos- 
mic rays. 

Measurements of the cosmic-ray antiproton spectrum 
started in the late 1970s and now extend up to ener- 
gies around 30 GeV [l|, 0, [!, Q • The observations can be 
extremely well described by adopting a production mech- 
anism as indicated above and then propagating the an- 
tiprotons through the galactic halo by means of a simple 
diffusion model @ . Remarkably, the parameters of this 
model are essentially fixed by the spectra of other cosmic 
ray species, in particular the boron over carbon (B/C) 
ratio [||. With the satellite-borne PAMELA experiment 
0| being in orbit since June 2006 and already taking 
data, and the planned installation of AMS-02 [|| on the 
international space station by 2007, the near future will 
witness very accurate and detailed measurements of the 
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antiproton spectrum, in the case of AMS-02 up to en- 
ergies of around 1 TeV. It is therefore of great interest 
to extend previous analyses of the background spectrum 
to higher energies since this will provide an important 
test of the underlying diffusion model and thus of our 
understanding of the propagation of cosmic rays through 
the galaxy. Having this in mind, the first part of the 
present work is devoted to a careful discussion of the un- 
certainties in the expected background spectrum at high 
energies and, as we shall see, how the already available 
data help to constrain it considerably. 

Not only is a thorough understanding of the cosmic 
ray antiproton background an interesting issue in itself, 
it is also a mandatory prerequisite for any attempt to 
spot possible exotic contributions. The main reason to 
expect these is the by now well-established existence of 
dark matter (DM) which, from the most recent measure- 
ments of the cosmic microwave background, contributes 
about 23 % to the total energy content of the universe 
Q. While the nature of dark matter is still an open 
question, weakly interacting massive particles (WIMPs) 
that arise in many extensions to the standard model of 
particle physics (SM) are very plausible candidates; ther- 
mally produced in the early universe, they automatically 
acquire the right relic density today As has first 

been noticed in [TTJ] , the self-annihilation of these parti- 
cles in the galactic halo could then leave an imprint in 
the antiproton spectrum. Traditionally, one has mainly 
focused on a possible contribution at rather low energies 
since first measurements seemed to indicate an excess 
of antiprotons below the peak at around 1 GeV 
With the improved statistics of follow-up experiments 
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and a better understanding of the production mecha- 
nisms, however, the evidence for such an excess disap- 
peared [HI]. In this article, we will therefore rather focus 
on possible exotic signatures at the high energies that 
soon will be accessible for the first time. To this end, we 
will compare the situation for three benchmark models 
of realistic dark matter candidates with masses in the 
TcV range; the lightest Kaluza-Klein particle (LKP) in 
models with universal extra dimension [Til . [l5| and the 
lightest supersymmetric particle (LSP) in the limit where 
it is a neutralino with a very high Higgsino or Wino 
[ItJ fraction, respectively. For earlier work that explic- 
itly studies the effect of dark matter annihilations on the 
high-energy antiproton spectrum, see e.g. [l8L fl9l. I20L l2ll] . 

This article is organized as follows. In Section [Til we 
briefly review the two-zone diffusion model that is used to 
propagate charged particles through the galactic halo; as 
mentioned before, this model provides an excellent fit to 
the low-energy data. We then continue in Scction[TTTJwith 
a detailed discussion, both numerically and analytically, 
of the expected antiproton background for energies up 
to 10 TeV. Section IIVI introduces the benchmark dark 
matter models that then will be used to discuss possible 
signatures for exotic physics in the antiproton flux and 
the prospects for their detection. In Section [V] finally, 
we give a short summary of our results and present our 
conclusions. 



II. PROPAGATION OF CHARGED PARTICLES 
THROUGH THE GALAXY 

Whatever the mechanism responsible for their pro- 
duction, charged cosmic rays subsequently propagate 
through the galactic magnetic field and bounce off its 
irregularities - the Alfven waves. The resulting parti- 
cle transport is well described by space diffusion with a 
coefficient 



K(E)=KqPK & 



(1) 



which increases as a power law with the rigidity 1Z of 
the particle. In addition, since these scattering centers 
move with a velocity V a ~ 20 to 100 km s _1 , a second 
order Fermi mechanism is responsible for some diffusive 
re-acceleration. Its coefficient Kee depends on the par- 
ticle velocity (3 and total energy E and is related to the 
space diffusion coefficient K(E) through 
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In the case of nuclei and antiprotons, one also has to 
take into account that a certain energy loss rate 6 loss (E) 
is induced by adiabatic and Coulomb energy losses, as 
well as by ionization. Finally, galactic convection wipes 
cosmic rays away from the disk with a velocity Vc ~ 5 
to 15 km s _1 . 



After these preliminaries, and equipped with all the 
necessary notation, we may now write the master equa- 
tion for the distribution function ip = dn/dE in space 
and in energy as: 

d z {Vc V>) ~ KAip + d E {b losB (E)ip - K EE {E)d E ^} = q . 

. (3) 

This equation applies to any species - protons, antipro- 
tons or positrons - as long as the rates for production q 
and energy loss b loss (E) are properly accounted for. This 
general framework - summarized in Eq. ((3)) - may be 
implemented within the semi-analytical two-zone model 
which is extensively discussed in 0, Q and whose salient 
features are briefly recalled here. 

In this approach, a steady state is assumed and the 
diffusive halo is pictured as a thick wheel which matches 
the circular structure of the Milky Way. The galactic 
disk of stars and gas - where primary cosmic rays are 
accelerated - lies in the middle. It extends radially 20 
kpc from the center and has a half-thickness h of 100 pc. 
Confinement layers where cosmic rays are trapped by dif- 
fusion lie above and beneath that disk. The intergalactic 
medium starts at the vertical boundaries z = ±L as well 
as beyond a radius of r = R = 20 kpc. Notice that the 
half-thickness L of the diffusive halo is not known and 
reasonable values range from 1 to 15 kpc. The diffu- 
sion coefficient K is the same everywhere whilst the con- 
vective velocity is exclusively vertical with component 
Vc(z) = Vcsign(z). The galactic wind - produced by 
the bulk of the disk stars like the Sun - drifts away from 
them along the vertical directions, hence the particular 
form assumed here for Vc ■ The spallation of the interstel- 
lar gas by cosmic rays takes place in the disk. It leads to 
the production of secondary species like boron or antipro- 
tons. The latter is a background against the neutralino 
annihilation signal and will be discussed in section HTT1 

The diffusive halo is axisymmctric and the cosmic ray 
density vanishes at the radius R — 20 kpc. This con- 
dition is naturally implemented by the following series 
expansion for ip: 



ip (r, z, E) = J2 p i (*> E ) J o («< r/R) 



(4) 



The Bessel function of zeroth order Jo vanishes at the 
points at. The radial dependence of ip is now taken into 
account by the set of its Bessel transforms Pi(z, E). The 
source term q may also be Bessel expanded into the corre- 
sponding functions E) and the master equation Q 
thus becomes 



d z (V c Pi)-Kd 2 z Pi+K 



R 2 



P, 



(5) 



+2hS(z) d E {b loss (E) P t - K EE {E) d E P t } = Q % (z, E) . 

Here, energy loss and diffusive reacceleration arc confined 
inside the galactic disk - which is considered infinitely 
thin, hence the presence of an effective term 2hS(z). 
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The form of the source term Qi(z, E) that appears in 
Eq. ([5]) depends on the type of the particle species that 
is considered. In the case of antiprotons, the following 
mechanisms can in principle contribute: 

• Antiprotons may collide elastically on interstellar H 
and He. Because they are preferentially scattered 
forward, however, such interactions are innocuous 
and will be disregarded. 

• Antiprotons may also annihilate on interstellar H 
and He. This leads to a negative source term 
— Tp ln ip, where the annihilation rate r™ is de- 
fined as 

pann ^ann „ i -.arm „ /£>\ 

: p = CT pH w p n H + o-p Hc Vpn H o- (o) 
The annihilation cross section cr^ 11 " is borrowed 

p ri 

from [25|, [2^] and we have multiplied it by a fac- 
tor of 4 2 / 3 ~ 2.5, taking into account the higher 
geometric cross section, to get 0p5a e - The aver- 
age hydrogen and helium densities in the galactic 
disk have been set equal to tih = 0.9 cm -3 and 
tihc = 0.1 cm~ 3 , respectively. 

• The annihilation of DM candidate particles 
throughout the Milky Way halo generates pri- 
mary antiprotons. The corresponding source term 
qp lm (r, z, E) will be discussed in section [TV] No- 
tice that annihilations take place all over the diffu- 
sive halo. 



• The latter is not the case neither for secondary an- 
tiprotons - which are produced as high energy pri- 
mary nuclei impinge on the atoms of the interstellar 
medium inside the galactic disk - nor for tertiary 
antiprotons which result from the inelastic and non- 
annihilating interactions which these particles may 
undergo with the same atoms. Antiprotons may 
actually collide on a proton at rest and transfer 
enough energy to excite it as a A resonance. This 
mechanism redistributes antiprotons towards lower 
energies and flattens their spectrum [l3j . 

The rate for the production of secondary antiprotons, 
qp c (r,Ep), is discussed in more detail in Section HTT1 for 
tertiary antiprotons, it is given by: 

qf(r,Ep) = 

£ (4 - %) »h 4 i>p(r, EL) dEp 

-o-pH^px (Epj n H vp ^p(r,Ep) , (7) 

where the inelastic and non-annihilating differential 
cross-section in this expression can be approximated by 

Here, Tp is the initial antiproton kinetic energy. In order 
to take into account elastic scattering on helium, one sim- 
ply has to replace the hydrogen density by tih + 4 2 / 3 % e . 



With all these source terms specified, the full expression for the master equation describing the (Bessel transformed) 
antiproton distribution function Pi(z,E) becomes: 

d z (V c Pi) - Kd 2 z P % + K Pi + 2h5(z)d E {b loss (E)P t - K EE (E)d E Pi} = 

- 2 h 8{z) Tf" Pi + Qp™ (z, E) + 2 h S(z) {Qg. + } . (9) 

Integrating this relation along the vertical axis z - in particular through the infinitely thin disk - finally leads to a 
diffusion equation in energy which the Bessel transforms Pi (0, E) fulfill: 

AiPi(0,E) + 2hd E {b loss (E)P t (0,E) - K EE (E)d E Pi(0,E)} = 



2h {Q^ + Qfi} + 2^ dzQ^(z,E) e~ 2K ^(z) . (10) 



I 

The coefficients A4 that appear in the above expression tions Ti(z) are defined as 
are given by 



ME) = V C + 2hTf n (E) + K(E) 5 4 coth , 



(L-z)\/ S mh{^-L 



(11) In this work, we solve Eq. (fTT))) by following the method 
where Sf — (Vc/K) 2 + (2 a,/i?) 2 , and the vertical func- explained in the appendix B of 
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A completely different approach to describe the an- 
tiproton propagation through the diffusive halo relies on 
the existence of a Green function Gp. Such a function 
translates the probability for an antiproton produced at 
point S(xs-,ys, z s) to travel to the observer located at 
point M(x, y, z). The antiproton energy spectrum is 
given by the convolution of the Green function Gp with 
the production rate qp 



1>(M,E) = Jd 3 x s G p (M^S,E)qp(S,E) 



(13) 



Energy loss, diffusive reacceleration and tertiary produc- 
tion are inefficient above a few GeV. For our purposes, we 
may therefore neglect these processes when deriving the 
antiproton propagator Gp. The Milky Way is now pic- 
tured as an infinite slab of half-thickness L with a gaseous 
disk in the middle at z = 0. The antiproton propagation 
is invariant under a translation along the horizontal axis 
x or y. The master equation ([3]) needs still to be solved 
along the vertical direction z with the condition that Gp 
vanishes at the boundaries z = ±L. The construction of 
the Green function for antiprotons is inspired from the 
positron case - see in particular section 3 of [22| - with 
the difference that the antiproton energy does not change 
and that time is integrated out. Explicit expressions for 
Gp may be found in the appendix C.3.a of [23j. Because 
we are interested in the flux at the earth - i.e. within the 
disk, at z = - the integral (fT3|) may be recast into 

i>(xe = r Q ,y Q = 0,z Q = 0,E) = 

4vr / dz s / r s dr s Gp(© <— S, E) qp(S, E) ,(14) 
Jo Jo 

where the antiproton propagator simplifies in that case 
into the scries 



Gp(© <-£,£) = 

^ 1 



e -z s /r... +°° 



27rK{E) ^ C n M0)Mzs)Ko\j;* 



.(15) 



Here, the vertical functions <p n are given by 

4> n {z) = sin{e„(l- J)} , (16) 
where the coefficients are solutions to the equation 

in = nir - tan _1 (p^„) . (17) 

At a given antiproton energy E, two specific scales 
may be singled out; the scattering length r^ 1 = 
hTp in (E)/K(E) and a convective scale r" 1 = 
V C /2K(E). The parameter 



1 

P 



L 



L 



(18) 



then serves to compare these two. At high energies, above 
~ 100 GeV, r s and much larger than L due to 



the greatly enhanced diffusion coefficient K(E). As a 
consequence, we have p> 1 and thus as (n — 1 /2) it. 
The scale G„, defined by 



C'n _ ^ 1 / sin^„ 



(19) 



tends to L in this regime where p is large, i.e. when diffu- 
sion takes over disk annihilations and galactic convection. 
In the argument of the modified Bcssel functions of the 
second kind K in Eq. (fTS")) . the ratio r/L is multiplied 
by a factor JT^, where 



= e 2 



L 



(20) 



which reduces to in the pure diffusive regime at high 
energy. 

As demonstrated later in more detail, the Bessel ex- 
pansion and the propagator approach give similar results. 
The small discrepancies between the two methods will be 
discussed in section lTvl but we may already now safely use 
Eq. ()15[) to get a deeper insight into antiproton propa- 
gation. In particular, we see that at high energies, the 
diffusion coefficient K(E) in front of the expansion (fT5|) 
is the only term that depends on the energy E. In the 
following Section, we will return to this observation when 
we discuss the expected antiproton background at high 
energies. 



III. THE ANTIPROTON BACKGROUND 

Secondary antiprotons are produced by the spallation 
of the interstellar gas of the Milky Way disk - mostly hy- 
drogen and helium - by impinging cosmic ray primaries. 
In the case of the interactions between cosmic ray protons 
and hydrogen atoms, the source term takes the following 
form: 



Q s p ec (r,E p ) 
~d~K 



(21) 



|s p —> £p j n H v p ip p (r, Ep) dE p 



The proton and helium cosmic ray fluxes at the Earth 
have been measured accurately and have been borrowed 
from Q. For each cosmic ray model, they have been 
consistently propagated backward in the whole diffusive 
halo in order to yield ip p (r, E p ) and ipH e {r, E nuc = E p ). 

For kinematic reasons, the production rate peaks 
around a few GeV; this is because the proton energy must 
be larger than a threshold of E® = 7m p . In the galac- 
tic frame, the differential production cross section that 
enters in the previous relation is given by the integral 



da 
dK 



2 7rfcf 



E 



d 3 a 



d(- 



(22) 
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where 9 denotes the angle between the momenta of the 
incoming proton and the produced antiproton. In the 
center of mass frame, which drifts with a velocity /3cm = 

1 /2 

{(.Bp — m p )/(E p + m p )} ' with respect to the galactic 
frame, the antiproton energy cannot exceed a value of 



E- 

p,max 



9 m" 



in: 



1/2 

where = {2m p (E p +m p )} is the total energy of 
the reaction. In Eq. (f2"2"| . the energies E p and Ep have 
been fixed and the angular integral runs from 9 = up 



(23) 



Et 



cos 9 n 



0CM kr 



p,max 



r 



CM 



(24) 



The Lorentz invariant differential cross section 
Ep (d 3 a / d 3 kp) depends on the antiproton rapid- 
ity y = tanh -1 (fcp \\/Ep) and transverse mass 
?ti^ = my + kp ± ; it has been parameterized ac- 
cording to [25L [2r| . In order to avoid numerical problems 
in the computation of (|22|) , it is convenient to change the 
integration variable from # max to mx for small angles 

^raax • 

So far, we have discussed the spallation of hydrogen by 
high-energy protons. Let us now turn to helium, which 
is also present in both the cosmic radiation as well as 
in the interstellar material of the galactic disk. As re- 
gards the spallation of interstellar helium by cosmic ray 
protons, the derivation of the differential cross section 
da/dEp(pHe — > p) follows the same lines as explained 
above. Because the incoming proton interacts with a 
single nucleon of the target nucleus, the center of mass 
frame of the collision is the same as before. The Lorentz 
invariant differential cross section - which is the only 
new ingredient - has been parameterized in [2?| as a 
function of t/s, kf\ , and the ratio xr = E$l/E$l 

v 5 Pj- 1 - n p/ p,max 

The total inelastic cross section for p + A collisions is 
taken from [28[; a boost is then necessary to cope with 
the collisions of high-energy helium nuclei on interstel- 
lar hydrogen. An antiproton with rapidity y in the 
reference system of the hydrogen target has a rapidity 
—y' = coshT 1 (E nuc /m p ) — y in the reference frame of 
the impinging helium nucleus. The integral (|2ip runs now 
over the energy per nucleon E nnc of the incident cosmic 
ray He. Finally, because helium amounts to a fraction 
of ~ 0.1 with respect to hydrogen, we expect He + He 
interactions to contribute only a few percent to the sec- 
ondary antiproton flux. Inspired by |29j, we have multi- 
plied the differential cross section da /dEp(-pHe — > p) by 
the effective number rj of nucleons inside the impinging 
He nucleus that participate in the reaction. It is given 
by the ratio of the average number < n mt > of nucleon- 
nucleon interactions over the average number < nx > of 
target nucleons involved. Both < m n t > and < tit > 
depend on the center of mass energy. At y/s —10 GeV 
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SECONDARY SPECTRUM 
Solar Minimum with F = 500 MV 



10 100 
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FIG. 1: The various contributions to secondary antiprotons 
from the spallation of the interstellar medium by cosmic rays. 
Here, we took the 'medium' configuration of propagation pa- 
rameters from Table [I] and Tp OA denotes the antiproton ki- 
netic energy as measured at the top of the atmosphere. For 
reference, we also show the existing low-energy data on the 
antiproton flux at the top of the atmosphere 0, 0, 0, 0j- 



for instance, we have < 7ij nt >= 3.68 and < tit >= 2.64, 
thus leading to n = 1.39. This value increases up to 
n = 5.12/3.25 = 1.58 at yfs = 1 TeV, not too far from 
the simple scaling factor of A 1 ! 3 — 1.59 which we have 
used in our code. 

In Fig. [TJ the various contributions to the secondary 
antiproton flux, $ = vip/Air , from the spallation of in- 
terstellar H and He by cosmic ray protons and alpha par- 
ticles are presented together with the existing low-energy 
data. Galactic propagation parameters corresponds to 
the 'medium' configuration of Tabic [I] to be discussed 
later. Solar modulation has been implemented through 
the force field approximation [30l ]. with a Fisk poten- 
tial c/)f of 500 MV corresponding to the minimum of 
solar activity during which the observations have been 
performed. Above an energy of a few tens of GeV, the 
antiproton flux decreases like a power law. If the scal- 
ing violations of the differential production cross sec- 
tion da/dEp were negligible, the secondary source term 
<7p c (r, E) would have the same energy dependence as the 

impinging cosmic ray flux $ = vip/Air oc E" 1 . As dis- 
cussed at the end of section|TTl the antiproton propagator 
Gp scales at high energy like 1/K oc E~ s . We expect the 
antiproton flux to have a typical spectral behaviour like 
Gp x $ oc E^ 1 ^ 6 . The cosmic ray proton and helium 
fluxes have been borrowed from [j| where a fit of the 
BESS [1H and AMS [H| data is featured. The spectral 
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FIG. 2: Theoretical uncertainties in the secondary flux in antiprotons, taking into account the whole range of propagation 
parameters that is allowed by the existing B/C data, again featured together with the existing low-energy data. In the right 
panel of this figure, we have plotted Tp&p in order to better illustrate the expected near Tp 3 scaling of the flux at high energies. 
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8 


K [kpc 2 /Myr] 


L [kpc] 


Vc [km/s] 


V a [km/s] 


max 


0.46 


0.0765 


15 


5 


117.6 


med 


0.70 


0.0112 


4 


12 


52.9 


min 


0.85 


0.0016 


1 


13.5 


22.4 



TABLE I: Typical combinations of diffusion parameters that 
are compatible with the B/C analysis @; as demonstrated 
later, these parameter sets correspond to minimal, medium 
and maximal primary antiproton fluxes, respectively. 



index 7 is found to be equal to 2.72 for protons and to 
2.74 for helium. Would hadronic interactions be scale 
invariant, the secondary antiproton flux of Fig. [T]- for 
which 5 = 0.7 - would drop like ~ E~ 3A . The actual 
spectrum is slightly harder with an E~ 3 3 energy depen- 
dence. As a cross-check, we also took for comparison the 
parameterization of [27j (instead of (2f| [26[) for the an- 
tiproton production cross section through pH-t pX. In 
this case we find qualitatively the same spectrum, with 
a slightly harder spectral index of 3.2 instead of 3.3. 
This corresponds to an uncertainty in the background 
flux which is much smaller than the one induced by the 
not fully determined propagation parameters - see the 
discussion below. 

The semi-analytic treatment of cosmic ray propaga- 
tion that is based on the Bessel expansion ^ is a con- 
venient framework to derive the theoretical uncertainties 
associated to the various parameters at stake - namely 
Kq, S, V a , Vc and the thickness L. The space of these 



propagation parameters has been extensively scanned [6j 
in order to select the allowed regions where the predic- 
tions on B/C - a typical cosmic ray secondary to pri- 
mary ratio - match the observations. Several hundreds 
of different propagation models have survived that crucial 
test. The propagation parameters are thus only loosely 
constrained by the cosmic ray nuclei abundances so far 
observed. The same conclusion has been reached inde- 
pendently in Ref. [33j ] with the help of a fully numerical 
code [HI in which the convective wind Vc increases lin- 
early with vertical height z. However, the B/C ratio 
could not be accounted for when both galactic convec- 
tion and diffusive reacceleration were implemented at the 
same time, a problem which our Bessel treatment does 
not encounter. 

The yellow band shown in Fig. [2] is the envelope of the 
secondary antiproton spectra computed with the set of 
~ 1, 600 different propagation models found in [|| to pass 
the B/C test. This band comprises the theoretical uncer- 
tainty in the determination of the secondary antiproton 
flux. It is confined by the 'minimal' and 'maximal' con- 
figurations of Table |TJ As a first observation, we notice 
how narrow the uncertainty strip is between ~ 10 and 
100 GeV. The PAMELA and AMS-02 collaborations will 
thus be able to highlight even small spectral deviations 
in that energy range. Above ~ 100 GeV, the yellow band 
widens as a result of the energy dependence of the dif- 
fusion coefficient K: from the B/C analysis, the spectral 
index S may take any value between 0.46 and 0.85; its 
spread A<5 = 0.4 thus translates into a factor of 10 8 ~ 6 
of uncertainty on the antiproton flux at 10 TeV - two 
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decades above the energy where the yellow strip is still 
the thinnest. In fact, we see this expectation confirmed 
in the right panel of Fig. [2j where we can read off a 
ratio of 1 to 6 (at 10 TeV) between the minimal and 
maximal antiproton flux expectations. This large uncer- 
tainty in the secondary antiproton background at TeV 
energies may look depressing. One should keep in mind, 
however, that PAMELA and AMS-02 will considerably 
improve the measurements of the cosmic ray nuclei abun- 
dances with a determination of the B/C ratio to a better 
accuracy and over a wider energy range than available 
so far. This will translate into improved constraints on 
the propagation parameters and eventually into a thin- 
ner uncertainty strip in the panels of Fig. [2] The an- 
tiproton spectrum itself will also be measured up to a 
few TeV in the case of AMS-02. Once it is compared 
to the cosmic ray proton and helium fluxes, the spectral 
index S should be better determined. Finally, we ex- 
pect the LHC to improve the accuracy of the antiproton 
production cross sections of the various nucleus- nucleus 
interactions at stake. 

When going to very high energies, a word of caution is 
nevertheless mandatory. The antiproton flux rapidly de- 
creases with energy and becomes vanishingly small above 
1 TeV. Even with the AMS-02 large acceptance, the sta- 
tistical error on a measurement of the antiproton back- 
ground at ~ 800 GeV is comparable to the present theo- 
retical uncertainty as featured in the left panel of Fig. [7] 
the vertical error bars of the highest energy data point 
extend over the yellow uncertainty strip. In fact, above 
10 TeV we expect not more than 0(1) antiprotons per 
m 2 and year, posing a formidable (if surmountable) chal- 
lenge even for follow-up experiments. In this article, we 
therefore restrict our analysis throughout to an energy 
range < 10 TeV. 
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NFW 97 [37] 
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5.38 


21.75 


Moore 04 [38J 


1 
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1.16 


2.54 


32.62 


Moore 99 [391 


1.5 


3 


1.5 


1.06 


34.52 



TABLE II: Parameters in Eq. (J26J) for the halo models that 
we consider here. Scale radius r s and density p s are strongly 
correlated with the virial mass of the galaxy [4t| and we adopt 
the values obtained in [4l[ for the Milky Way. In the case 
of the Moore 99 profile, DM self-annihilations set an upper 
bound to the maximal possible density and we follow the usual 
prescription of imposing a cutoff radius inside which pcDM 
is assumed to be constant [i^]. When the DM distribution 
is cuspy - for 7 > 1 - we smooth it so as to keep the total 
number of annihilations constant - see the text for further 
details. 



The source function can thus be disentangled into an 
astrophysical part, given by the dark matter distribution 
Pcdm, and a particle physics part that comprises the 
remaining information. The first one is usually subject 
to much greater uncertainties than the latter and will be 
dealt with in the following subsection. In section IIV B[ 
we will then address the particle physics part by dis- 
cussing three benchmark scenarios with DM masses in 
the TeV range, with their respective expressions for the 
annihilation cross section and branching ratios. The frag- 
mentation functions AN* /AT, finally, can be obtained 
from Monte Carlo programs like Pythia [35[. We use 
here the tabulated fragmentation functions of the Dark- 
SUSY package (36|, which are based on a large number 
of Pythia runs (10 6 per annihilation channel and mass). 



IV. ANTIPROTONS FROM DARK MATTER 
ANNIHILATIONS IN THE GALACTIC HALO 



The source function q£ rim (T, r) describes the number 
of primary antiprotons per unit time, energy and volume 
element that are produced with a kinetic energy T at a 
given position r in the galactic halo: 



prim 



(T,v) 



1, 

2 Wain 



Pcdm(v)\ 2 TB fdNf 



f 



AT 



Here, (cr a 



(25) 

is the DM self-annihilation rate, pcdm the 
DM halo density and m the DM particle's mass. The sum 
runs over all annihilation channels /, where B* are the 
branching ratios and AN* /AT the fragmentation func- 
tions into antiprotons, respectively. The factor of 1/2 
has to be included whenever the DM particle is its own 
antiparticle and therefore only annihilates in pairs (which 
is the case for the DM candidates that we will consider 
below) . 



A. The distribution of dark matter in the Milky 
Way 

1, Halo profiles 

The Milky Way halo distribution is only poorly con- 
strained by direct observations and its form must there- 
fore be inferred from N-body simulations of gravitational 
clustering. Due to the restricted resolution of these simu- 
lations, however, the innermost slope of the density pro- 
file can only be obtained by an extrapolation of the be- 
haviour at larger radii (i.e. r > 0.1 kpc) and is thus 
bound to a considerable amount of uncertainty; further- 
more, it is still an (numerically) unsolved issue how to 
correctly include the effect of baryons during the gravita- 
tional collapse. To account for this situation, we will here 
for illustrative purposes restrict ourselves to a choice of 
four different profiles that basically span the whole range 
of reasonable halo models with respect to indirect dark 
matter detection prospects; the isothermal sphere, the 
Navarro- Frenk- White (NFW) H3 and the Moore profiles 
[38l . l39l ]. Like most of their relatives, these halo models 
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can be parameterized as 
PCDM(r) = A 



r 



(26) 



where the corresponding parameters are summarized in 
Table [IIJ Note that at the sun's distance to the galactic 
center, r = Rq = 8.5 kpc, one recovers for all profiles a 
local halo density of po = Pcdm(-Ro) ~ 0.3 GeV cm -3 . 

Our semi-analytic treatment of cosmic ray galactic 
propagation relies on Bessel expansions like (|3]). The 
function Jo (ct!j r/ii) probes details as small as A ~ 
27r Rj on where on is the i-th zero of the function Jo and 
R denotes the radius of the diffusive halo. In the case of 
an NFW or Moore profile, the correct description of the 
central cusp - where the annihilation rate diverges like 
r -27 - would imply the necessity of an infinite series of 
such Bessel functions. This would lead to unacceptable 
CPU time and to a lack of numerical accuracy. Following 
the method explained in [43| . we have renormalized the 
DM central distribution without modifying the absolute 
number of DM particle annihilations. Actually, within a 
small sphere of radius r c , we have replaced the density 
profile (|26p which diverges like p(r)/p c = {r c /r) 7 where 
p c = p(r c ), by the milder distribution 



p(r < r c ) 

Pc 



= 1 



2 7r , , o / 7T r , , 

_ (e _l) sin 2(_)^ 



where sin c (a;) = sin(x)/x. That renormalized density 
leads to the same number of DM annihilations as the 
actual cusp provided that 



3 - 27 ' 



(28) 



in the case of NFW 97 and Moore 04. Note that for the 
Moore 99 cusp, the actual annihilation rate diverges log- 
arithmically and we have to impose a cutoff on the DM 



density, p C DM 



< 



Pa 



This upper bound arises from 



the increasing DM self-annihilation rate in regions of en- 
hanced DM densities and is given by: 



m 



(29) 



where an age of t = 12 Gy has been assumed for the 
central DM concentration. This cutoff translates into a 
factor of 



i 



Mooro 99 



= 2 In 



Pc 



(30) 



In practice, we have set r c = 500 pc and pushed the 
Bessel expansion up to the 200-th order. We have accel- 
erated its convergence by using the method explained 
in [44|. Because the antiproton propagator Gp that 
connects the solar system to the galactic center varies 
smoothly over that region, a smaller core radius r c would 
not appreciably change our results. 



2. Boost factors 

Following the paradigm of hierarchical structure for- 
mation, it is very likely that the DM distribution exhibits 
substructures ("clumps"), i.e. small inhomogencitics su- 
perimposed on the smooth background profiles described 
above. On very small scales, e.g., one expects a cutoff in 
the cold dark matter power spectrum [45j and a consid- 
erable number of the first gravitationally bound objects, 
with a mass around this cutoff, might have survived until 
today [46] . Another interesting proposal is the existence 
of a population of intermediate mass black holes (IMBHs) 
in our galaxy, each of them being surrounded by a DM 
mini-halo with a rather steep profile, i.e. a "mini-cusp" 
[47l ]. The effect of halo substructures is in any case to 
enhance the flux of antiprotons - or any other primarily 
produced particles - at earth, as compared to the case of 
a smooth DM matter distribution. Such an enhancement 
immediately follows from the fact that (p 2 ) is generally 
larger than (p) 2 , with the difference increasing with the 
degree of inhomogeneity, and is often encoded in a uni- 
versal boost factor b. However, as has been stressed in 
[22| , this boost factor is actually not so universal in that 
it - given a DM distribution pcdm ~ depends on both the 
particle species and the energy range under consideration 
(see, e.g., [2(| for a similar observation). 

As an example, and as a practical application of the 
formalism developed in [22| . let us now determine the 
boost factor for antiprotons in the "type B" IMBH sce- 
nario of [47I ] - which is representative of a class of mod- 
els in which black holes originate from massive objects 
formed directly during the collapse of primordial gas in 
early- forming halos. A complete discussion of the an- 
tiproton signal and of its statistical properties is not in 
the focus of this paper and the reader is rather referred 
to the thorough analysis of Ref. [24| for details. As a 
sneak preview of that work, we derive here the effective 
boost factor 



B r 



< *p , IMBH > 



(31) 



p , smooth 



for an average IMBH population of "type B" . A smooth 
halo density pcdm yields a flux 



$p, — cc J d^ s G p (M^S,E) {^^y 



whereas the IMBH contribution is given by 



(32) 



p,IMBH 



>cx Jd 3 x s G p (M^S,E)£ s 



s P (r) • 
(33) 

In Ref. [24|, the average number n sp of mini-spikes per 
unit volume has been fitted with the results of ~ 200 dif- 
ferent Monte-Carlo realizations - each containing eventu- 
ally a hundred objects - obtained in Ref. [13] by evolving 
an initial distribution of IMBHs orbiting in the galactic 
halo and by taking into account close encounters and 
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FIG. 3: Boost factor for an average Milky Way halo population of IMBHs as a function of the antiproton kinetic energy. The 
two panels on the left do not include the contribution from tertiaries, energy losses and diffusive reacceleration. The long 
dashed lines have been obtained by integrating directly the antiproton Green function Gp over the diffusive halo. In the right 
panels, tertiaries, energy losses and diffusive reacceleration have been taken into account. An isothermal sphere - top panels - 
as well as an NFW profile - bottom panels - have been considered for the DM smooth distribution. 



their associated tidal disruptions. The enhancement of 
the DM annihilation rate inside each mini-spike is de- 
scribed by the source boost factor 

U = f (^W) 2 dr , (34) 

J DM cloud I P0 J 

where po is some value of reference generally set equal 



to Pcdm(-Ro)- In a region of size R sp around the IMBH, 
the DM density is described by a power law r -7 / 3 . The 
annihilation rate diverges at the center and we impose a 
cutoff on the DM density pcdm Pann ~ see Eq. (p?9")) 
and Ref. 0- The effective volume £ sp depends on 
the DM particle mass m and self-annihilation cross sec- 
tion (cr anll u). We have used here an average spike ra- 
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dius of R sp = 2.84 pc and an average surface density of 
<$Pcdm(-Rs P ) = 48.51 M Q pc -3 - see once again Ref. [24| 
for details. 

In the upper panels of Fig. [3l the isothermal sphere of 
Table |TT] corresponds to the smooth distribution of ref- 
erence with respect to which the enhancement of the 
antiproton signal is computed. A DM particle mass 
m = 1 TeV and a self-annihilation cross section (cr ann w) = 
3 x 10~ 24 cm 3 s _1 lead to a mini-spike effective volume 
of £sp = 1.2 x 10 5 kpc 3 . In the lower panels where the 
isothermal sphere is replaced by an NFW profile, a DM 
particle mass of 1 TeV and a self-annihilation cross sec- 
tion of 3 x 10~ 26 cm 3 s _1 have been assumed, hence a 
value of £ sp = 3.1 x 10 6 kpc 3 . In the left panels, cosmic 
ray propagation does not include the contribution from 
tcrtiaries, nor energy losses and diffusive reacceleration. 
Switching off these processes allows to perform the calcu- 
lation of B e g by the two methods presented in section HT1 
The convolution over the diffusive halo of the antipro- 
ton propagator Gp with the square {pcdm(i")/po} 2 °f 
the smooth DM density profile (J2BJ) and with the IMBH 
source term £, sp n sp (r) leads to the long-dashed curves. 
Alternatively, we have directly included the IMBH source 
term ^s P n sp (r) in our Bessel code to generate the solid 
lines. The injection spectrum of antiprotons cancels out 
in the calculation. The agreement between these two dif- 
ferent approaches is astonishingly good. The discrepan- 
cies do not exceed a few percent in the worst case. This 
occurs in particular for 'maximal' propagation where the 
range of the antiproton propagation is large and allows 
to probe the radial boundaries of the diffusive halo. In 
the Bessel approach, the cosmic ray density vanishes at 
r = R. This condition is not implemented in the prop- 
agator approach even if the contribution of the outward 
regions is not taken into account. This leads to an overes- 
timation of B Q g as is particularly clear in the upper-left 
pannel. At first glance, the boost factor depends very 
weakly on the propagation model. As the antiproton en- 
ergy increases, so does the propagation range. In the 
case of the isothermal sphere, the population of galactic 
IMBHs - whose number density is given on average by 
n sp (r) - tends to contribute more to the signal than the 
smooth distribution. As the propagation range increases, 
more and more mini-spikes come into play. This results 
into the increase of B ff in the upper-left pannel for the 
'medium' and 'maximal' propagation configurations. On 
the contrary, in the case of an NFW profile, as soon as 
the propagation range reaches the galactic center and its 
strong DM annihilation site, the boost factor drops as is 
particularly clear for the 'medium' configuration of the 
lower-left pannel. In the right panels, tertiary produc- 
tion, energy losses and diffusive reacceleration have been 
switched on with the consequence of smoothing the varia- 
tions of the boost factor with energy. In Fig. [31 the values 
of £ sp have been adjusted in order to get a boost of ~ 10 4 

though the DM particle mass and self-annihilation cross 
section are reasonable. Because the distribution of mini- 
spikes inside the galactic halo is not unique, the boost 



factor B e g is subject to large variations - typically a fac- 
tor of a few exceeding 10 at low antiproton energies - de- 
pending on the actual IMBH distribution. For a full anal- 
ysis of the associated boost uncertainty we refer to [24| 
- but would like to mention already here that we expect 
this variance to decrease significantly at the high energies 
considered here. 



B. Dark matter masses at the TeV scale 

Searching for possible primary contributions to the an- 
tiproton spectrum at very high energies, we have to fo- 
cus on the annihilation of rather heavy (i.e. at least TeV 
scale) DM particles. While such high masses can appear 
in a variety of models, we will restrict ourselves in the fol- 
lowing to the most popular and most commonly studied 
scenarios, i.e. supersymmetry and extra dimensions. In 
this section, we introduce the DM candidates that arise 
in these scenarios and in each case briefly describe the 
main features that are relevant in our context. First, 
however, note that the usual WIMP relation, 

a 3-10- 27 cm 3 s- 1 

"WIMP" ~ ; : , (.oO) 

(av) 

in general rather prefers slightly smaller masses. To bal- 
ance the reduced annihilation cross sections one expects 
for higher masses, one has thus to rely on the presence 
of effective coannihilation channels with other particles 
during the freeze-out process. In accordance with this 
observation, the mass splitting between the DM particle 
and the next-to-lightest non-standard model particles is 
generically very small in all the examples described be- 
low. 

Let us now start with the case of supersymmetry, 
where the lightest stable supersymmetric particle (LSP) 
provides an excellent dark matter candidate (Toj . In most 
models, it is given by the (lightest) ncutralino, which is a 
linear combination of the superpartners of the gauge and 
Higgs fields, 

X = X°i = NuB + N 12 W 3 + N 13 H° + N 14 H° . (36) 

While the neutralino is often a gaugino, with a large Bino 
fraction and a mass of a couple of hundred GeV or less, 
the hyperbolic branch/focus point region of minimal su- 
pergravity (mSUGRA) typically exhibits very heavy neu- 
tralinos with a large Higgsino fraction [l6j . From the re- 
quirement that Higgsinos should give the right relic den- 
sity, their mass has to be around 1 TeV [48[. For these 
high masses, the neutralino is an almost pure (anti-) sym- 
metric combination of the two neutral Higgsino states - 
in which case the annihilation cross section into Z (as 
well as into Higgs boson) pairs vanishes exactly and that 
into quarks is usually heavily suppressed by multi-TcV 
squark masses in the propagator. For these reasons the 
annihilation into W bosons typically clearly dominates; 
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depending on the actual neutralino composition, how- 
ever, a significant fraction can also go into Z bosons or 
heavy quark pairs. 

Another interesting situation arises in the case of a 
neutralino that is almost a pure Wino, as expected for 
example in anomal y m ediated supersymmetry breaking 
(AMSB) scenarios [49|]. For Winos, the preferred mass 
from relic density requirements is peaked at about 1.7 
TeV [HI (a pure Wino state obtains the right relic density 
for m ~ 2.2 TeV; taking into account non-perturbative 
effects, this mass is further increased by about 600 GeV 
[5(il]). Non-perturbative, binding energy effects then re- 
sult in greatly enhanced annihilation cross-sections to- 
day, when the neutralinos have very small galactic ve- 
locities [13] • In this limit, heavy Winos annihilate al- 
most exclusively into gauge bosons (as a side remark, the 
annihilation into photons is also significantly enhanced 
w.r.t. the non-perturbative result, leading to promising 
prospects for an indirect detection in terms of gamma 

rays 0MH3). 

The third example of a TeV scale dark matter can- 
didate we want to consider here is that of the lightest 
Kaluza-Klcin particle (LKP) in models with universal 
extra dimensions (UED) [lj], where all standard model 
fields are allowed to propagate in a higher-dimensional 
bulk. After compactification of the internal space, these 
additional degrees of freedom appear as towers of new, 
heavy states in the effective four-dimensional theory; the 
stability of the lightest of these states, i.e. the LKP, is 
guaranteed by the existence of an internal Z2 symmetry 
(called KK parity) that derives from higher-dimensional 
translational invariancc. Taking into account radiative 
corrections to the KK masses, the LKP is expected to be 
well approximated by the , the first KK excitation of 
the weak hypercharge boson [15|. Detailed relic density 
calculations show that it can account for the required 
dark matter density if the compactification scale (and 
thus the mass) lies in the range 0.6 < m B (i) < 1.4 
TeV [53[ , mainly depending on the standard model Higgs 
mass; deviations from the minimal scheme for calculat- 
ing the radiative mass spectrum weaken the upper bound 
on the compactification scale to about 2 or 3 TeV [54| . 
The main annihilation channels of the B^ that are rele- 
vant for antiproton production are those into quark pairs 
(about 35 % in total); the annihilation into gauge and 
Higgs bosons is of the order of 1 % each and thus sub- 
dominant. 

Finally, let us stress once more that in all these exam- 
ples of heavy WIMP candidates, TeV masses appear in a 
very natural way. Furthermore, one may easily find even 
higher masses for e.g. neutralino dark matter as soon 
as one leaves minimal prescriptions for supersymmetry 
breaking [55| . However, since the annihilation cross sec- 
tion in general scales as (er a nnv) oc m~ 2 , the number of 
primary antiprotons is suppressed as 

prim —4 /ot\ 

<7p oc m . (37) 
This means that even if the background flux scales 



roughly as T~ 3 (see Section UTT]) . it is very unlikely that 
one will be able to discriminate a primary antiproton 
contribution from DM particles with masses significantly 
higher than 1 TeV - unless one allows for artificially high 
boost factors [6(| . Since also the total rates become very 
low far beyond 1 TeV, we do not consider such high DM 
masses in the following. 



C. Future detectability 

After having introduced typical scenarios with high 
DM masses, let us now choose as benchmark models a 
1 TeV pure Higgsino, a 1.7 TeV pure Wino and a B^> 
with corresponding masses (having in mind, of course, a 
direct comparison between the LKP and the LSP). In Ta- 
ble IIII1 we summarize the main properties of these DM 
candidates. Before we now proceed to study possible 
imprints in the antiproton spectrum for our benchmark 
models, we note that the cases of a pure Wino [2l[ and 
the LKP [2(| 5H have already attracted some attention 
before. 

With the source function specified in Eq. (|2"B"|) , we can 
easily apply the formalism described in Section|TT]to com- 
pute the corresponding primary flux of antiprotons. Let 
us start by showing in Fig. [4] and Fig. [5] the results for 
our benchmark models, for various halo profiles and dif- 
fusion parameters. Note first that the dependence on the 
halo profile is as usual much less severe than in the case 
of e.g. gamma rays, where the flux usually differs by sev- 
eral orders of magnitude for the profiles that we consider 
here; the simple reason for this, of course, is that the vast 
majority of antiprotons do not reach far enough through 
the diffusive halo to probe the very inner region of the 
galaxy, where the differences between the halo profiles 
are most pronounced. The dependence on the diffusion 
parameters, on the contrary, is rather strong - although 
the B/C constraints basically fix the secondary antipro- 
ton flux. This effect has been described before [5j,l56| and 
can be attributed to the fact that primary and secondary 
antiprotons mainly probe different regions of the halo. It 
clearly illustrates the need for cosmic ray data that are 
both more accurate and span a larger energy range, each 
of which would greatly increase the predictability for pri- 
mary contributions to the antiproton flux. 

Unfortunately, the expected primary components are 
usually smaller than the background fluxes in secondary 
antiprotons even for favourable diffusion parameters - 
with the striking exception of the 1.7 TeV Wino, where 
the resonantly enhanced annihilation cross section may 
allow for a spectacular signal in the range of a few 100 
GeV already for very conservative assumptions about the 
DM distribution (i.e. no boost factor at all). We can 
thus confirm the claim of [2l[ that Wino DM exhibits 
very promising observational prospects in terms of pri- 
mary contributions to the antiproton flux. 

However, let us now recall from our discussion in the 
proceeding Section that we actually do expect primary 
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TABLE III: Our benchmark models for studying possible primary contributions to the antiproton spectrum at high energies; 
m is the DM particle's mass (in TeV), (a^nnv) its annihilation rate (in 10~ 26 cm 3 s~ 1 ) and the remaining columns give the 
branching ratios into the annihilation channels relevant for p production (in percent). The corresponding values are typical for 
high Higgsino (LSP1.0) and Wino (LSP1.7) fractions of the neutralino; for the latter we take the non-perturbative expressions 
from [1T]]> while for the former we have calculated the annihilation cross section of a pure (anti-)symmetric Higgsino into W 
bosons and choose to neglect other annihilation channels (see text for further details). In the case of the LKP, the quoted 
values for the branching ratios, as well as (cr ann ?;)m 2 , are actually very insensitive to the parameters of the model (see, e.g., 

0)- 



fluxes that are boosted with respect to what is shown in 
Figs |4] and [5] As we have emphasized, even rather high 
boost factors are feasible and appear in several realistic 
scenarios. To illustrate this point, we plot in Fig. [6] the 
primary fluxes for our benchmark models - as before for 
an NFW profile and the diffusion parameter configura- 
tions of Table [H but this time in the presence of IMBHs. 
All of a sudden, the situation has drastically changed and 
a plethora of new features in the antiproton spectrum 
seems to lurk just behind the currently reachable ener- 
gies ! In fact, with primary fluxes that much enhanced 
one even has to start worrying about compatibility with 
the existing low-energy data; such an analysis, however, 
would be beyond the scope of the present work and war- 
rants a dedicated future study. 

Finally, let us address the question whether the next 
generation of experiments will be able to distinguish be- 
tween the models that we have presented here. To this 
end, we consider in Fig. [7] for each model a boost fac- 
tor that normalizes the maximal deviation from the sec- 
ondary flux to that of the 1.7 TeV Wino case. For the 
high energies that we are interested in, the precision of 
cosmic ray flux measurements is essentially limited by 
statistics (see also the remark at the end of Section IIIII) . 
For comparison, we have therefore included in Fig. [7| 
the statistical error after 3 years of data sampling by 
PAMELA and AMS, respectively, provided that these 
experiments would measure an antiproton spectrum as 
induced by LSP annihilations. As we have seen before, a 
1.7 TeV LSP can easily be distinguished from the back- 
ground (even when taking into account the full uncer- 
tainty in the spectrum of secondary antiprotons) already 
by PAMELA. In fact, this is true without having to in- 
voke any boost factors at all. In order to discriminate the 
spectra of the other benchmark DM candidates at a sim- 
ilar confidence level, one would need boost factor from 
about 150 to 500. While such boost factors may seem 
rather high from a traditional point of view, we have 
stressed above that, e.g., in scenarios with DM mini- 
spikes around IMBHs one would naturally expect even 
higher values. 

When it comes to the actual discrimination between 



different DM candidates, however, we realize that the 
prospects are less promising; given the current uncer- 
tainty in the secondary flux, as well as the expected sta- 
tistical errors in the data, neither PAMELA nor AMS 
will be able to distinguish between different types of an- 
nihilating WIMPs (i.e. LSP vs. LKP). A determination 
of the WIMP mass, on the other hand, will be possi- 
ble to a certain extent - at least when a clear drop in 
the spectrum becomes visible (for the DM models con- 
sidered here this could be the case once the AMS data 
are available). Note also that there appears in the spec- 
trum a certain degeneracy between the WIMP type and 
its mass, putting a principle limit on the accuracy of any 
possible mass determination: for a given mass, the LSP 
produces an annihilation spectrum that is very similar 
to that of an LKP, apart from being slightly shifted to 
higher energies. An LKP with a somewhat enhanced 
mass (by about 10%) would therefore feature a spectrum 
that is almost indistinguishable from that of the LSP. 



V. SUMMARY AND CONCLUSIONS 

Upcoming experiments like PAMELA and AMS will 
extend the upper range of accessible antiproton cosmic 
ray energies considerably, from hitherto around 30 GeV 
to well above 100 GeV. In this article, we have considered 
the galactic antiproton spectrum at high energies and 
performed a detailed analysis of what these - and suc- 
cessive - experiments are expected to see. To begin with, 
we have presented the first discussion of the non-trivial 
behaviour of the antiproton background at high energies. 
Our analysis has been performed in the framework of a 
two-zone diffusion model, whose parameters can be ex- 
tracted from the spectra of other cosmic ray species (the 
B/C ratio, in particular) in such a way as to give an ex- 
tremely satisfying explanation of the observed low-energy 
antiproton data (notice, however, that this is not the con- 
clusion reached in [5^ where galactic diffusion is assumed 
to be anisotropic). We have shown that the diffusion pa- 
rameters obtained in this way determine, in fact, almost 
uniquely the spectrum of secondary antiprotons up to 
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FIG. 4: The primary flux in antiprotons, as compared to the background in secondaries, for the set of benchmark models 
specified in Table IIIII From bottom to top, the different curves correspond to the isothermal sphere, NFW, Moore 04 and 
Moore 99 profiles, respectively; for the diffusion parameters we adopt the 'medium' configuration of Table U 



energies of around 100 GeV. For larger energies, the un- 
certainty in the expected spectrum slowly increases, up 
to a factor of about 6 at 10 TeV. The main reason for 
this uncertainty lies in the fact that the B/C analysis 
does not sufficiently constrain the spectral index S of the 
diffusion coefficient; this uncertainty will be considerably 
reduced once experiments like PAMELA and AMS-02 im- 
prove the quality and range of existing B/C data. Based 
on both numerical as well as semi-analytical results, we 



furthermore predict a simple power-law scaling for the 
high-energy secondary spectrum, above slightly less than 
100 GeV, with a spectral index in the range 3.1 to 3.5. 

In the second part of this paper, we have investigated 
the possibility of a contribution from primary antipro- 
tons to the high-energy part of the spectrum, originat- 
ing from the annihilation of DM particles in the galactic 
halo. To this end, we introduced a set of benchmark 
models with high DM masses, motivated by DM candi- 
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FIG. 5: Same as Fig. [4] now for an NFW halo profile; the diffusion parameters are varied from the 'minimal' to the 'maximal' 
configurations of Table [I] 



dates naturally arising in theories with supersymmetric 
or extra-dimensional extensions to the standard model 
of particle physics. We have shown that the expected 
primary antiproton signal takes in all of these cases a 
spectral form that is sufficiently different from the back- 
ground to clearly discriminate it against the latter if the 
corresponding total fluxes are high enough. For this lat- 
ter requirement to be satisfied, one usually has to invoke 
boost factors of 0(100) that at first sight seem rather 
high (with the interesting exception of a Wino-like DM 



candidate, see also [21(). As a side- result of the present 
work, however, we have also derived the expected boost 
factors in the IMBH scenario of [13] and demonstrated 
that in this case, in fact, they are even higher than what 
is needed for a clear attribution of the observed signal to 
a DM (as opposed to secondary antiproton) origin. We, 
furthermore, note that the most recent, high-resolution 
simulations of gravitational clustering [53] hint at a con- 
siderable DM fraction to be distributed in clumps of 
highly enhanced DM densities, which potentially would 
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FIG. 6: Same as Fig. [5] now with an NFW profile and a typical population of IMBHs in the galactic halo. 



lead to large boost factors, too. 

When it comes to the possibility of distinguishing be- 
tween different DM candidates, however, the antiproton 
spectrum does not turn out to be very well suited and 
other indirect DM detection methods (gamma-rays, in 
particular) are probably more promising in that respect. 
The mass of the annihilating DM particle, on the other 
hand, may eventually be determined, to a fair accuracy, 
from the antiproton spectrum alone; a comparison with 
the result of other DM searches would then give an impor- 
tant, independent piece of information. It is also worth- 



while to study in this context the cross-correlation of an- 
tiproton signals with DM-induccd features in the spec- 
trum of other charged cosmic ray species (for the LKP, 
e.g., this would be a pronounced peak in the positron 
spectrum that necessarily has to appear if the antipro- 
ton spectrum shows an LKP-induced distortion [2fj, [58| 
- a feature that is not shared by LSP DM candidates). 

To conclude, measurements of the galactic antiproton 
spectrum are an important testbed for our understand- 
ing of the propagation of charged particles through the 
galaxy and the diffusive structure of the Milky Way halo. 
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FIG. 7: The yellow band shows the expected antiproton background for the full range of allowed diffusion parameters. In 
the left (right) panel, we have added the primary and total fluxes for our 1 TeV (1.7 TeV) benchmark DM models of Table 
Mil assuming an NFW halo profile and a 'medium' set of diffusion parameters; for illustration, and to better compare these 
models, we have adopted a boost factor of 2 (1050, 330, 270) for the case of an LSP1.7 (LKP1.7, LSP1.0, LKP1.0) dark matter 
candidate. For the LSP, we also include the expected statistical error after 3 years of data sampling by PAMELA and AMS-02, 
respectively. 



As a means of indirect DM detection, antiprotons are 
probably not the single-most promising species; we would 
like to stress, however, that one should, rather, focus on 
the complementarity of different approaches, to which 
antiproton measurements would then have the potential 
of providing an important contribution. 
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